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ABSTRACT 


Herein is presented a derivation and computational 
formulation for a modified Kolomogorov-Smirnov test. This 
test extends the hypothesis testing and confidence limits 
advantages of the Kolomogorov-Smirnov test to data which 
is censored beyond a predetermined number of observations. 

A Vwstange of the computer program nsedmiaaceae calenus 
lations of significance levels and the resulting signifi- 
cance levels for specified parameter values for the modified 


test are included. 
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I. INTRODUCTION 


im statistacalwapplicatwons, it is frequently desired 
to determine, with some specified confidence, whether a 
yandom sample is fromea hypothesized distribution. Sim- 
itm; satteissfrequently desirable to establish confidence 
limits on the distribution function using sample observa- 
tiems. Mere are numérous procedures which accomplish 
these objectives when the sample is uncensored. The 
Kolomogorov-Smirnov test is applicable to both of these 
objectives and has the advantages of being distribution 
free amd simple to apply. 

In situations wherein the data is censored beyond a 
pReaetcamumned number of faiiumes, the complete sampie is 
not available since the values associated with sample ob- 
Se€rvations greater than the censoring point are unobser- 
wable. In such a situation, One Still mikoihtecesire to 
conduct a “Kollomogoroyv-Smirmnoy type'’ test over only the 
ordered uncensored sample observations. Such a modified 
Kolomogorov-Smirnov test is derived herein. To the know- 
ledge of the author, this test has not been previously 
available. 

The modified Kolomogorov-Smirnov test is proved to 
be distribution free. Furthermore, a computational for- 
mulation is presented for calculating significance levels 
and the significance levels are tabulated as a function 


Of Sample size and critical value. 





Proofs and derivations contained herein do not speci- 


fically deal with the usual problems with events of prob- 


ability zero, however they can readily be modified to 


include such events. All computer work performed in 


connection with this paper was conducted at the United 


States Naval Postgraduate School on an IBM 360 computer. 


The following listing constitutes the description 


and/or definition of notation utilized herein: 


Notation 


S 
n 


ene = SX 
aes 


>) - 


Descripiion 


Cumulative distribution function of 
a random sample of size n. 


Hypothesized distribution fume F Lome 


Number of observations over which the 
modified test 1S run. 


Level or test. 

Sample size. 

Cr remeal values. 

Critical value associated with a test 
run over the first f ordered observation 


from sample size n. 


Ordered observed random sample from dis- 
tm DUtTOn F . 


2 (get | st) = EG) |< De ) 


X<X ¢ 





II. KOLOMOGOROV-SMIRNOV ES | 


Kolomogorov established that the distribution of the 
Sita Stic 


supmiesmed) - FCO | 
xX 


is independent of the parent population distribution F if 
F is continuous, and he derived its limiting distribution. 
Massey [3] and Miller [4] present methods by which signif- 
icance level a can be computed for a given critical value 


d such that . 
PiG@ Sue | cae3) - F(x) | > d) = a. 
Xx 


This relation between significance level versus critical 


La act + + ~~ 
the basis for applying the 


thi 


value and sample size 1 
Kolomogorov-Smirnov test. Further information on the der- 
ivation, tabulation and use of the test 1s contained in 
References [1], [2], [3] and [4]. 

The essence of the Kolomogorov-Smirnov test 1S ClieEs 
for a sample distribution function (S_) and a continuous 


parent population distribution function (F), the sStatlorme 


sup |S. - F | 


n 
is independent of the distribution of the parent population. 
Thus, given that the population has distribution F, the 
probability that the maximum absolute difference between 
F and the sample distribution will exceed a specified 


critical value is dependent only upon the sample size and 





the critical value specified. The test itself then re) Sa ul 
ple to conduct since it entails only the determination of 
this maximum difference and whether it is less than the 
specified critical value D.- Furthermore, as stated by 
Kendall and Stuart (1, pg 457) "we may reverse the pro- 
cedure of testing for fit and use dD, to set confidence 
limits for a (continuous) distribution function as a 
whole." 

The Kolomogorov-Smirnov test is thus a simple method 
for either testing goodness-of-fit or establishing con- 
fidence limits for a distribution function from a sample 
when a continuous parent population distribution is hy- 


pothesized. 





III. EXISTENCE OF THE MODIFIED KOLOMOGOROV-SMIRNOV_ TEST 


The Kolomogorov-Smirnov goodness of fit test establishes 
a band Si + d such that the probability is (1-9) that eine 
true distribution function F lies entirely within this band, 


foe}. iat is, 


B( supe | SOs; hie.) | > d ) =a; 
x 


where ao is the probability that the test rejects the hypoth- 
esis that F is the distribution of X, given F is the distri- 
bution of X and where d is the critical value for testwsize 
a. Let this test be denoted by T. Then T rejects iff 

ae | S(x) > F(x) | > d for some x. 

The modified Kolomogorov-Smirnoyv test. @uoposed hereun 
is to apply the Kolomogorov-Smirnov test using d* for eras 
ical value and applying the test only to the first £ orden 
statistics to yield a significance level oa. iat. 1S); 


P( SUP | Seo) = Hoe) > ole) eee 
X<X ¢ 


Let this test be denoted by T*. 
Let A denote the event that rejection occurs for some 


Xs i < f when using test T, given that test T rejects. Then, 


with each test using critical value d*, 


P(T* rejects) P(T nejyects ir (teweeurotlag sain firstmeesorden 


statistics/rejection occurs) 


= P(T rejects) P(A) < P(7T rejects) = a. 





However, P(T rejects) is a continuous monotone de- 


creasing function of d. Therefore, for each d there is a 


d* (d*<d) such that 
P(T rejects using d) = P(T rejects io iimemacn b Un 


P(T* rejects using d*) 


which establishes the existence of the test Ra 





IV. PROOF THAT THE MODIFIED KOLOMOGOROV-SMIRNOV TEST 15 
DISTRIBUTION FREE 


Lec Yi» ¥ 


continuous distribution F. Define 


go-to ty be an ordered random sample from a 


w= sUP | S(y) - FO) | 
ver: 


whe re niece) = inebe US? Yy Veen et 


X, = F(¥;) which may be written "X, (Y)" 


and 


W* = SUP | S*(x) - x | 
XX =: 


Then Xy> X, 7-7 ok, is (almost surely) an ordered sample from 
the uniform distribution over the irmenvar (0,1). 


However, 


W 


SUP | ae) - F(y) | 
re 


sup } | max {ki F(Y,) < F(y)} - 2 FO) | 
ye 


: SUP 1 | max {k: F(¥,) < F(y)} - 2 Fy) | 
{y:F(y)<F(¥,)} 


= SUP 1 | Max ik: Xy See > 2 X(y) | 
fy :X(y) <X¢} 


ou | max {k: X, < x} = mx | 
n Kom 
X$X ¢ 


= SUP | S*(x) - x | 
n 
x<X > 


_ ws 


10 





Thus W = W* and it follows that for any critical value D 
Pee eee PW > DI 


Since W* has a distribution not depending upon F, it follows 


that the test is distribution free. 


alg 





VY. DERIVATION OF CRITICAL VALUES (D, ew) 
9 


a 


Mie significance level, Pent? associated with a 
givem Critical value, De in? can be calculated by the fol- 
lowing procedure. This procedure was motivated by, #@and ais 
an analytical generalization of the procedure presented by 
Massey [3] to calculate Kolomogorov-Smirnov critical values. 

Since it has been shown in Section IV that Pe 6D) is 
independent of F provided only that F is continuous, it 1s 


sufficient to consider only the case 


BIO) = 


as 
be 
Fy 
2 
A 
yas 
1A 
an 


j--2 
He 
rin 
ra 
Vv 
ae) 


Divide the interval (0,1) into n parts by the points 
1/n, 2/n,---,(n-1l)/n. The step ssDU Mie SL Se rises by jumps 
of exactly 1/n. Thus, in order to be inside the band 
Bie) + k/nm at x = i/n, S_(i/n) must be one of the lattice 
points whose ordinates are (i-k+1)/n, (i-k+2)/n,---,(itk-1)/n. 


Let a, be an n dimensional vector of non-negative intewe rs 


“77 40 . ) 


a. = (GQ,.., GG, oe Ge : 
1 ( ieee 2G aeapse 


suemetiat 


fA 


j= 


Suppose that the step function 2 stays inside the band by 


means of Oe 4 of the observations from F falling within tine 
? 


az 





iomcrere) (( ja) f/m, j/n), 7 = 1, 2,---,n. Under H,, the 


probability of this happening is given by the multinomial 





law as 
tz, 
n! ia? J> 

we ee ee es Ta 
ee 1. gure @ G's to i! On 4 {= 

= il a! 

ar oi! mn 

ial J,1 


as the probability of the step function staying within the 


bemd F + se” 1s given by 


= 1 ny .. me i 
a? -) n . >. n 
i 


; i Th i ° i 
A 5 Oe. ite ; ‘oar 1° 
= J» ae! J> 


where the summation is over all i for which a, is associated 
with She Within the band. There exists a one-to-one corre- 
spondence between ta; } and the paths of Sat 

Define a = (a5, Gor ~ 71), where the elements of a 
mepresent all paths of Sn: Bac clement sot a Lrepresemims 
ewener a path for which S09 remains entirely within the 
band tor X < X¢ or a path for which SO) does not Temadimecn— 
Simemy within the band for x < Xp- 

Let I = {O55 S109 remains entirely within the band 


fon xX < Xe}. Then by the multinomial law 


en il 
P ( SUP | S09 =) X. | < ae = Fa y > ‘ 
ae a a iel a. .! 


j-iaeee 


Consider all possible values of M1 k (that is Se ae (ee 


nm). “aen.ecither 


re 





(a) malik ieee 


OT 


(b) 1k ee 


and one of the following mutually exclusive cases must hold 
(refer to Figure 1). 
Ate. if A, 2 f and S60) is entirely within the 
bamd for x < Xe, then o,e€1 regardless of the 
values of on he --- i The contr@bucion you 


a, to ee 1s 


Tis ! O ! ! 


Spas ee” ok 





nh 
j=2 aiek 


where the sum is over all a. k Suen ot dae 


1 Ga “ee 

faa. =| ae a, , 2 f and S09) 1ssmoet entimelyawathin thie 
band for x < x,, then a,.#l regardless of the value 
eae 9k? Toe Oy and therefore Oy does not con- 
tribute to the sum which composes Peas 

sc] J ae ls A. f and S09) is entirely within the band 
fom x < xX , then whether a,eI will depend upon 

a al | k 

Oo? ee k Aide tieweconcriputdon of a, to 
Pen?) 1s dependent upon Oo 1? Os ---, Ok" 

fee, at oe f and S09 is not entirely within the 
band for x < x, , then a, gI regardless of >? 


1 Ok 
Oe ke? 77 Gt i and therefore ty does not contri- 


(De 


bute to the sum which composes ey 
2 
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Consider next a, = Ny a ie ---, ok? such that 


ay , =? Bll. Then either 


Or 


Ce ee 2ak 2 FE 


(b) a < f 


lp # 72,8 


and one of the following mutually exclusive cases must hold 


Brefer to Figure 1). 


AZ . 


AZZ. 


Bl2. 


i “1k 4 Oo kk = & “and S029 1s entirely within 
the: band’ fer x < Xe, then a,el regardless of @ehe 
values of ae O4 ke? ies On k and the contribu- 


tron of Oy to oa 1s 


where the sum is over all a. sSlteh that .2 > (qe 
yk j= as 


-— ek A 
ant 1k + Oo k > £ and SOx) iS NOt entimely, waun- 
im the bend for xe< Xe» then ap el regamdless "08 
Ok? = SA On k and Oy does not contribute to the 
sum which composes Pe 76D) 

1¢ 1 Ok + 19k < £ 4am S40) is entirely within 
thie band: femex < x 


Tak 22k 
will depend upon Ok? An ke "772 WH and the con- 


2 , then whether a,el 


tripwet10on of a, to Pe 7 6D) is dependent upon 


OS ae? ’ oa ae’ 


ee 





pee. if O41 k +h x < £ and SE (C3), is not entirely 


within the band for x < X, then 


+ @ 4 
1,k 2,k 
ap el and Oty does not contribute to the sum which 


composes an) 
Consider next o, * (04 4? hy 27728 Kk? such tivat 
+ E> Bll and =e Continue this procedure until 
the contribution from each Oss eo 1 ae P 76) is obtained. 
Then 

(1) oe = the sum over all ieI of the conta bis 


tions to Dee from Ge 








i 
ew. 
/ il = Oy .t---a, et 
— 
> at Lag 
i Ea 
" Tees o3.j | ae 
il 
= oe 
0. 1 j‘°  snee j 
n-2 jopAl.n-2 2 b 
1 mn : 
n 1 Ag O-1 jap Al net 
ee ee. 
ae 
ee ne 


n 
where the inner sums are over all a. -. such that ole 
1,) ; Lo 1,) 
Let Ukr) =) goes over a such that ‘ a. = Yr. 
a, }---O,: i=l ea 


Then U(k,r) satisfies the difference equation 


16 





, 1 : 
(2) U(m+1,3) = a Gay Ue) 50,1, 


with initial conditions of 








UCOei iy = ai 
and 
O(0 4) = 20 eae CO - 
Thus,.by substituting (2) into feb); Py 7 6D) is of the 
b 
form 
P @) - ee L Dine bs nea 
eT n a : Med 
n i) 
oy pall 
\- : uC 
s = =s - U(n-2, n-a, ,-a, ;) 
i oe 13 N20) 
oj 
Oy 5 a2! 
n= 2,3 
U(2,n-a; 5A j ee a i Ae na §e 
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Upper limit of band 


A21 Al2 as) Al4 


A21 Alz Al3 Al4 


: a 
4/ 





f/n=3/9 A21 Aulg2 Al3 Al4 
lower limit 
" eA — of band 


£/ 5 





0 1/5 275 NS 5 1 
Parameters: 
n= 5 
D = 275 
f= 3. 
Pmeure 1. Classification of Lattice Points. 
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VI. COMPUTATIONAL FORMULATION FOR CALCULATION OF SIGNIFI- 
CAMCE LEVELS 
The foregoing derivation can be utilized to calculate 
the significance level associated with a given critical 
value for the modified Kolomogorov-Smirnov test. However, 
the procedure would be cumbersome and not amenable to com- 
putation. A more tractible computational extension is dis- 
cussed in what follows. 
Massey [3] has derived the set of 2k-1 difference eque- 


tions 


. 


U det) U, (a) , = 1, -- "ee 


Gstl-h)! Dose 


with 


U, (a) =O ith>mt tk 


U; (0) = 0 Gee i = & 
U, (9) =i for b= & 
which yields 
E/¥a yp = 78 
Pack/vm ) = FY Uy(n) 


for the probability that the sample distribution function 
Le.) stays entirely within the band F + k/n. 

Referring to the (0,1) interval divided into n equal 
parts and indexed with l through n, the I oe valves repre- 
sent values assigned to the lattice points within the band 
F + k/n (refer to Figure 2). 

The index m refers to the ordinate index f the lat- 


tice point. Thus, the values for the lattice points with 


19 





VI. COMPUTATIONAL FORMULATION FOR CALCULATION OF SIGNIFI- 

CANCE LEVELS 

The foregoing derivation can be utilized to calculate 
the significance level associated with a given critical 
value for the modified Kolomogorov-Smirnov test. However, 
the procedure would be cumbersome and not amenable to com- 
putation. A more tractible computational extension is dis- 
cussed in what follows. 

Massey [3] has derived the set of 2k-1 difference equa- 


tions 


hazel . 
Il _ — ' 
U; (m+1) yy Gti-h)! U,, (m) | lace BASSI 
i= 1 


with 
U, (m) = Oren oes fs 
U; (0) = 0 for i #k 
U, (9) = | @or i = k 


which yields 


P_(k/vn ) Hatt, (on) 


for the probability that the sample distribution fumeeron 
ty) stays entirely within the band F + k/n. 

Referring to the (0,1) interval divided into n equal 
partsmeand Indexed with 1 tnmoien Nemec Oe (us, values repre- 
sent values assigned to the lattice points within the band 
B® k/m (refer to Figure 2). 

The index m refers to the ordinate index »f the lat- 


tice point. Thus, the values for the lattice points with 


Js, 





ordinate (m+l) depend only upon the values of lattice points 
within the band with ordinate m and upon the critical value 
(k/n) and sample size (n). 

U, (n) then represents the value associated with paths 
On Sa which remain entirely within the band F + k/n and pass 
DimoucheAattice point (N,n).. Consequently U, (n) + (n!/n") is 
the probability that Sh will remain entirely within the band 
|e hoe 

Massey's approach can be modified to provide similar 
difference equations for the calculation of the significance 
Eemelwassociated with critical valwe (k/n), number of un- 
censored failures (f), and sample size (n) for the modified 
Kolomogorov-Smirnov test. The required modification is to 
peparce the réquirement thet 2 ey Temaan Sneirely wena 


F + k/n'' with the requirement that i remain entirely with- 


m F * k/n through the - 


Ordered Observation." 
The set of difference equations which results from this 


modification is 


j@ 
7 1 
ojgeert) = d (Cae igen 


Nol 
co 
i 
we 
=. 


where the upper limit on the sum is given by 


n-m+k-1 if £ < F + k/n 


j+l a f.> F + Wme(eeter tor Piper 


Additional computational efficiency can be obtained by 
recalling that the U; (m) values for each lattice point re- 


Sults from the oe components whick. determine specific paths 
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m 
for Sn: Thus, if .2, 4, 2 f then the first f ordered ob- 


i=1 
servations remain within the band F + k/n regardless of a,, 
i = m+l, ---, n. Consequently, the U; (m) values for these 
iieetce points can be relabled as Z. (m) and the difference 
equations summed only over the lattice points which repre- 
sent paths for which it is undetermined whether Sn remains 
within F + k/n for the first f ordered observations. Ref- 


erencing Figure 4, it may be seen that this yields the equa- 


tions 


a il 
U, (m+1) 2 Gale U, Gn) 


= 1 
2m 2) nec, GFI-R)T Uy Gm) 


where 5 represents the circled lattice points in Figure 4 
having abscissa values less than or equal to j. 
Consider the set of difference equations 
J 


OVUP (i+1,)) <a GeicayT OVUP SH) ; j=0 eee 
h=0 


with initial conditions 


OVUP(0,0) = 1 


OVUP(0,}) = 0 if j # 0. 


Orie) 1S an array which represents values associated 
with path segments for which Sn traverses i lattice points 
to the right while increasing in abscissa value by j lattice 


points. Consequently, the Zen values can be translated to 
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values at lattice point (n,n) by multiplying each Z 5 (m) 
value by OVUP(n-m, n-j-mt+f-1). The value of U, (2) is then 


the sum of these translated Z; (m) values. That 15s 


U, (n) - » Z; (m) OVUP(n-m, n-j-m+i-1) 
il pdb 72 


and 


P(F-k/n < Sn < F+ek/n) = U, (n) “1 


yields the significance level associated with the critical 
value k/n, k an integer, for the modified Kolomogorov- 
cmernev test. 

This computational procedure 1s specifically applicable 
to critical values of k/n, k an integer. However, the pro- 
cedure can readily be adapted to criticai values of k/ ime 
with k and m integers by dividing the (0,1) interval into 
m°n parts wate n parts. Moreover, the significance Weed 
associated with critical values not Specifically caleulared 
can be approximated by interpolation between critical values 
for which the significance levels were enllenlatved: 

This latter compuational formulation is the basis of 
the computer program presented as Appendix A. The signifi- 
cance levels as calculated from the program in Appendix A 
are presented as Appendix B for all combinations of the fol- 
lowing parameters: 

PSS 10, 15, 205. 2ommes0 
ee? 35 4 15 omar enmeo gtO 


(574 





Jesh eabnesy f4e 





Massey's Procedure. 


Zo 








7 
a 
7 


U-_ (2) Uaesp U.(4) 
5\@ 1% 3°"@ 


ri 
7 


wad 
Upeenue- Ueaan U., (4) 
anc A 3® é 


7 
7 





U,(2) U, (4) 
2s Loe 





—_ —_—_ —_—_— ar —_—_-  —_—_——l = -_— ——— 


7 
7 


phe OD 
ia 


Pesce. 5. 





U, (2) 


Modified Procedure. 
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U, (5) 














Z, (2) Z,(3) Z_(4) 
Sip oo 3 oo 









7 
Z4(2)_, 7233) 
oo =e 





Pe 
ATs 





Z., (3) 


N 


Figure 4. Modified Procedure Using Translation. 
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VII. VERIFICATION OF SIGNIFICANCE LEVEL COMPUTATIONS 


In order to verify the correctness of the computer pro- 
gram, a computer simulation was conducted. This simulation 
was based upon a population distribution which was uniform 
over the interval (0,1). The simulation was run for sam- 
ple sizes (n) of 5 and 10, with critical value (k/n) of 
6940, and 3000 replications. 

The computer program for the Simulation is listedaan 
Appendix C. The significance levels resulting from the 
simulations and the corresponding significance levels from 
Appendix B are shown in Table I. It appears that the simu- 


lated values and the calculated values are in agreement. 


re) 


This contention that the computed values are fapproximately) 
correct is further supported by comparing them with the 
values calculated by Massey [3]. These values are shown in 
ape 1. 

In view of the foregoing, it 1s concluded that the 
calculated significance levels (Appendix B) and the computer 


program from which they were calculated (Appendix Aj are 


accurate and contain no significant errors. 
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TABLE I 


SeeriCch | '  # OBSERVATIONS ! SIGNIFICANCE LEVEL 
VALUE : no ee : SIMULATED [CALCULATED 
40 | 5! 1 1.92330 1.92224 
; | 2 | . 8663 .87104 
3 | 3 "478670 ' 78752 
| | 4 | 47277 * 72960 
, | S ' 6923 ' — .69120 
40 7 10a 1 1.9937" 99395 
| | 2 9883 ' 698614 
. . 3 9803 ' .97906 
4 ' .9740 ' 97 59m 
' 
| S 6.9677, 96953 
7 | 6 1.9620 ' 96331 
. 7 1.9560" —.95596 
; . 8 r  .9463 | .94895 
| 9 1 69413 1.94365 
10 1.9387 ' Seeatoel 


a SS yy ai SS So eS i SS SSS SSS GS SSS SS ss 


a 





TABLE II 


SAMPLE ' CRITICAL ' ; 
SIZE (n) VALUE MASSEY APPENDIX B 
: ! 20 0384 03840 
40 6521* 69120 
60 : 9699 96992 
80 99935 999 36 
t t 
10 10 . 0004 00036 
20 eS 25128 
30 7331+ 72946 
40 9410 94101 
; 50 9922 99222 
460 9994 99943 
20 : 10 0238 02374 
: iE 2955 29553 
2 6473 64728 
j a: 8624 86237 
: 30 9569 95693 
: . 9892 98924 
: 40 9979 99787 
| 45 ) 9997 99967 
RE 1 20 7637 16367 
24 9057 90564 
1 28 9683 96832 
! a7 1 9911 99109 
2G 9979 99791 
t t 
30 . eed . nego 66289 
2000 } 8420 84202 
55 9359 93587 
2667 oy 74 97744 
3000 9931 99314 


t 


*significance difference exists between Massey's results and 
results within Appendix B. Further study has indicated that 
the results of Appendix B are correct. 
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VIII. APPLICATIONS OF MODIFIED KOLOMOGOROV-SMIRNOV TEST 


The modified Kolomogorov-Smirnov test is essentially 
a Kolomogorov-Smirnov test run over only the first f 
ordered observations. Therefore, the modified test 15 9ap- 
plicable to data for which the Kolomogorov-Smirnov test 
would be applicable except that the data is censored at 
some predetermined number of ordered observations or, for 
some other reason, it is undesirable or impossible to “een- 
sider datum beyond the specified .number of observations. 
When the modified test uses all n ordered observations, 
the modified Kolomogorov-Smirnov test 15s mentacal to tile 
Kolomogorov-Smirnov test. 

As with the Kolomogorov-Smirnov test, the modified 
test can be used as a goodness of fit test if the (contin- 
uous) distribution function (F) is completely hypothesized. 
Furthermore, just as the probability statement fommbhie 
Kolomogorov-Smirnov test can be inverted to yield a method 
of setting confidence limits (1, pg 457), the probability 
statement for the modified Kolomogorov-Smirnov test can te 
inverted to yield a method of setting confidence limits 
atadizing only the first f ordered sample observations. 
This attribute of the modified test can be useful In@ecn. 
sored life testing situations where data is censored at some 


predetermined number of ordered observations. 


z¢ 





IX. WORK REMAINING TO BE DONE 


Areas wherein it appears prudent to direct fUtieMte beet ils 
vestigation of the modified Kolomogorov-Smirnov test pre- 
sented herein are: 

(a) evaluation of the power of the modified test. 

(b) comparison of the performance of the modified 
test against other test procedures presented in the liter- 
ature. 

(c) restructure the test to provide critical value 
versus significance level value for censoring at a pre- 
determined point instead of a predetermined number ene allo 
servations. 


= 
aa 


° * a a oe ys a ns ~~ °° a 
It is conjectured that 1tem fc} can be accompiisnec 


ap) 


by utilizing the difference equations presented by Massey 
[3] except carrying the summation over lattice points Wren - 
in the band only to the censor point and then translating 
these values to the n/n lattice point by use of the 


OVER Cin) array presented in Section VI of this paper. 
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X. SUMMARY 


The Kolomogorov-Smirnov test 1s a well documented 
procedure for testing goodness of fit between a hypothe- 
sized distribution and an unknown distribution by uti- 
lizing a random sample drawn from the unknown distribution. 
The Kolomogorov-Smirnov test can be utilized to establish 
confidence limits for the hypothesized @istaibution funcs 
tion as a whole. The significance levels associated with 
given critical values are well tabulated and procedures 
for their calculations have been published in the litera- 
Etre . 


A derivation and computational formulation for a modi- 


bd 


fied Kolomogorov-Smirnov test has been presented. tL ex- 
tends the hypothesis testing and confidence limits advantages 
of the Kolomogorov-Smirnov test to one for censored data. 

The modified test is essentially a Kolomogorov-Smirnov test 
run only over the first f ordered sample observations. A 
listing of the computer program used in the calculations 

of significance levels of the modified test for given 
critical values is presented as Appendix A. The resulting 


significance levels for specified parameter values auc 


given in Appendix B. ; 


eal 





APPENDIX A: COMPUTER PROGRAM FOR CALCULATION OF SIGNIFICANCE 


= 
a 


jo 7 
arp} 


A 
= 


LEVELS 


LCULATIGN TF CENSOPED PRCBABILITLaS» CV=K/No K AN INTEGER 
DIMONSTAM Z(31,31), BVHE8(30,39), Peers) ) 
PATA OVUP/900*0.0/; FATT/1.0,30*%0.0/ 

MTAC = 31 

TOL. 'AC=2 Lo 

FACT(T’C) = FACT(IAC—-L) #FLOAT(IAC-1) 
10 COATING: 

CVUE(L,1) = 1-9 

AN 4Nn T= 25 

7 39 J=) 9 30 

Meee = 1+ 1 

aqvuce(tTs,J) = 0.0 

NO 290 NH=1,35 

muMe( ie. J) = IVE) + OVUP (1-1 .NH)/FSCT (J+1—-NH) 
BOe ON] TNT 
30 CONTINUE 
wore Corry tT MUE 

r= tN VUES (Is bk TS FOR ONER (1-1) .,AND UP (J-1l) 

SREEOR-, TO GO AVES T AND WP 3 CALL QVUPCIF1, J+1) 
ae 3200 N=5,329,5 
RIT: (6,5) ™ 

5 FOOMET(5X,"*N = ©, 13) 
JIN = N41 
Mad = “NO(19,N) 
DN 20 C K= ver 
WRITS (626) K 
& Smet T (5X, ** = *, 13) 
PO 110 M&=l,"! 
N00 45 T=1,31 
NT 44 1=1,31 
Z(T,5) = 0.2.9 
og Ce. Nit 
45 Cir) og iS 
1s 1) — ma 
MTA = MINO(M=+K, N41) 
me POO 1ee=2,- % 
MINGa = MAXO(1L. TA-K+1) 
We (1f..070. Nei) MINIGSN4+1 
MAXef = N+1 
D7 OC JASMIN MX 
7A tee, SO) = O.0 
NHMIM = MAXO(L, TA-¥) 
NHMVAX = MTINQ(K4+TC=—2 9 we @ fA) 
DQ 70 Miser TU Tomales 
Mofo | = Li tags) + Bil o= ee ey we 1 ( eo 
70 GaaipoMi'+ 
SO CUT hi: 

100 CONT IMIS 

ee OME CT GELeut Lien JF THE Z AND CVJIE VALUES ae 

PAUP TCDS NyKyNE WHERE M=SSAMPLE STZcy K/N=CRITICAL VALUE 

n NES# ORSEAVATIONS AVER WHICH THE TEST 1S RUN . 
TAMis} = MAXC(NE —-" #2 ,.1) 

TBMAX = MINQ(MF4K, ‘It1) 

Sum = 0.0 

aoe 50 Te=lP MIN, TPM XxX 

weeifh = N+) 

Tee (78 129. “+l ) AISMIN=NEL 

JBMAX = N+) 

Paet4o Jo= Jee, Bt 

Si = Si“ + 7( 1s JR) ZQVUP (N- 184+2, “N-J8+2) 

140 CONTINUE 

Toomey NUS 
Say = Sim = FAC ee ane = DAT (N) EN ) 

WRITS (6,690) ME, ey 

609 FARMAT(9X,13-' OSS -RVATIONS - - -— —',F9-5) 

Lie Basis 

200 We IMU 

S09) Coie = 
STI 
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APPENDIX B: CALCULATED SIGNIFICANCE LEVELS 


(n = sample size, f = number of observations, critical value = k/n) 


n= 5 
Et i 2 3 4 S 
f 6 7 g 9 10 
1 * emG7252 0.92224 0.98976 0.99967 0.99999 
: Be eee en ee a 
2 , 0.28000 0.87104 0.98496 0.99968 1.00000 
Meee ee a + 
3 ' 0.12160 0.78752 0.98496 0.99968 1.00000 
, Pe oe a 
4, 0.05760 0.72960 0.97824 0.99968 1.00000 
aa! ee ee. ee 
5 ' 0.03840 0.69120 0.96992 0.99936 1.00000 
iP SSS a. 
n = 10 
k 1 2 3 4 5 
£ : 6 7 8 Q 16 
1 ' 0.65132 0.89262 0.97175 0.99395 0.99902 
| 0.99989 0.99999 1.00000 1.00000 1.00000 
> 0.25320 0.81192 0.94152 0.98614 0.99771 
' 0.99977 0.99999 1.00000 1.00000 1.00000 
0 aire le w ee eee eee eee ee 
3 ' 9.09911 ~—dO. 68126 0.91515 0.97906 0.99668 
0.99972 0.99999 1 00000 1.00000 1.00000 
at eS eer oo) ee - = = aa eS hs 
i Olney 0.55947. 0.88112 0.97392. 9.99618 
0.99971 0.99999 1.00000 1.000001. 00000 
a ge - - = oa ms oe es 
PP amos 6d. 0.46168 0.83988 0.96953 0.9961 
0.99971 0.99999 1.00000 1 00000 - 00000 
ee ee a x - = = , fe _- =-« =e 
6 * 0.00635. 0.38556. 0.80282. 0.96331. 0.99596 
0.99971 0. 99999 1 - 00000 1.00000, 00000 
—. tae, = ae = = = - oo = eee 
7 eer 00265 0.32827 0177438 0.95596 0.99516 
j .99970 0.99999 I. 00000 1.00000 1.00000 
3 | 0.00115 0.28842 0.75368 0.94895 0.99391. 
0.99963 0. 99999 1.00000 1.00000 1.00000 
ge Wee 00054 0.26563 0.73821. 0.94368 0.99279 
0.99951 0.99998 1. 00000 1.00000 1.00000 
A. 5 ee Phe See = ee sie RS a Se = 
10 ' 0.00036 0.25128 0.72946 0.94101 0.99222. 
0.99943 0.99998 1.00000 1.000 °9 1.00000 


= 
= - - => = = = = -_ = = = = = = - -= -= oe -_= - - - = -_= = =~ = - = = - = - 


ao 





0 


0 


0. 


0 


0. 


“9 


0. 


0 


0. 


0.00015 
0. 
0. “00005. 
of 


“0. 00002. 
i 


“9. 


0. 00000 
Ue 


0.00000 


.64473 
0. 


aad 


uae 


00536 
QO. 


.00082 


0.88310 
99952 0. 


0.79514 
.99862 0 
09394 0.65514 

i. a QO. 


Ore 52460. 
.99640 v. 


01386 


0. 41916 
99549 UF 


Oe (33562. 
ae 9486 OF 


00208 
0. 


Oe 26985 
99450 On 


21822. 

Wr. 09413" QO. 
OF “17779 
99347 0. 


00032 
OF 


QO. 14631. 
— 0. 


QO. "12206 
ee gs 
Oy. 10388 
a 0. 
00001 QO. 09125 
0.98953 0. 


0. 03404 
mee he 
0. “07950 
2 


9888 0 


7 5 g 4 9 5 10 
0.96481 0.99045 0.99771 
99991 0.99998 0.99999 0.99999 
0.92578 0.97675 0.99379 
99975 0.99996 0.99999 0.99999 
0.88982 0.96212 0.98918 
99955 0.99993 0.99998 9. 99999 
0.84326 5 .caems,  Onaeeren 
999370. 999920. 99998 0.99999 
0.78520  O- “93326 0.98045 
99925 0.99991 0.99998 0.99999 
0.72783. 5 onmie 0.9768 
99919 0.99991 0.99998 0.99999 
0.67530 0.89308 0.97277 
99918 0.99991 0.99999 0.99999 
0.62879 5 ag270° 0 96810. 
99917 0.99991 0.99998 6.99999 
0.58861 n.gg45s 4—«..: 0.9633 
99914 9.999910. 99999 aa 99999. 
oe. ae 3946 0.95870 
999050. 999919. 99999 0.99999 
ee, “Op eee oman 438 
99888. 99990 0.99999 0. 99999 
0.50827. ~—«O. sie). oboumien 
99869 0.99988 0.99999 0.99999 
9.49415 0.81038 0.94777. 
99852 0.99986 0.99999 0.99999 
0.48377 1 poses 08 “94594 
99841 0.99984 9.99999 I. 00000 
0.47795 0.80275 0.94517 
99837 0.99983 0.99999 1.00000 
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re 2 3 4 5 
a! 6 7 8 9 10 
' + 
» * @.64151 0.87842 0.96123 0.98846 0.99681 
0.99919 0.99980 0.99994 0.99997 0.99998 
- © eta @« = - oS = aS = = - «- #« sees = Sa - 
2 0.24227 (0.78722 0.91799 0.97155 0.99112 
0.99751 0.99937 0.99985 0.99996 0.99998 
a = oe |) a ms oo = ~— Cabo oe = = 6 & SS -6) = = = eee 
se * 9.09163 0.64328 0.87783 0.95298 0.98398 
' 0.99520 0.99874 0.99971 0.99995 0.99997 
eee I = ey és = = - = ~— « « « =e 
4 0.03471 0.50956 _—O. 82612 (0.93455 0.97627 
Gno0253 MGs Ores 0.99990 0.99997 
= 2 = - oa = — - -©-= = = «= 
et 0.01317 0.40200 0.76163 0. 91448 0.96853 
0.98975 0.99718 0.99935 0.99987 0.99996 
5 | ee 5 = = - « «eS hee «= ) ae 
6 ' 0.00501 0.31713. 0.69734 0.88859 0.96080 
' 0.98708 0.99643 0.99919 0.99984 0.99996 
ea te See = a a oO = = a - = - =-= «= 
7 * 0.00191 0.25053. 0.63757 0.85935 0.95180 
0.98462 0.99577 0.99906 0.99983 0.99996 
g  ' 0.00073 0.19836 0.58334 5.82944, 0.94128 
0.98215 0. 99526 0.99898 = 0.99982 ).99996 
SO. ra er a = eee as a ee rah = - = ee 
9 ' 0.00028 0.15748 0.53468 0.80048 0.92996 
0.97940 0.99485 0.99894 0.99982 0.99996 
a Ee = Age = = ma — = - - 
10 ' 0.00011 0.12544 0.49130 (0.77337 0.91863 
0.97641 0.99443 0.99892 0.99982 0.99996 
= SE ee er ee es ‘a 2 = o = oo 
11 ' 0.00004 0.10032 0.45287 0.74859 0.90789 
' 0.97333 gxg9302 = «099888 0.99982. 99996 
eae) a” = - - 
12 ' 0.00002 0.08063 —_—O. "41906. «0.72646 OO. (89818 
0.97033 0.99329 0.99880 0. 99981 0.99996 
=e ae os en as a a . « «= 6 a= 
13 ' 0.00001 0.06519 0.38964 0. “70720 «0.88976 
' 0.96751 0.99255 q.s9g67 0.99980 0.99996 
ee oS) SS ee ee ime = = = - ©« e© - = 
14 ' 0.00000 0.05311 0.36443 0.69100 0.88268 
0.96488 0.99177 0.99850 0.99978 0.99996 
fs ee emt ae = = = — =« « © 8. 
15 ' 0.00000 _0. 94370 «0.54340 sO. .67795 0.87680 
' 0.96251 0.99101 0.99833 0.9987 99996 
2 at) 2 SS eee ee oe o = zk Sea = 
16 ' 0.00000 0.03646 0.32665 0.66783 0.87184 
! 0.96049 0.99036 0.99816 0.99972 0.99996 


ee, 
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kK! iE Z 3 4 5 


: 6 7 8 9 10 
17 ‘' 0.00000 0.03103 0.31434 0.65979 0.86785 
0, “ae 0. oF 0.99803 0.99970 0.99996 
rs ' O@. 00000 OF. 02726 0. 30558 0265 o54 0.864 
0. ge UL 0. eel 0.99794 0. 2208 at gees 
LO a0, 00000. 0. 02510 0. 29914 0. 64931 0.86312 
0.95716 0.98930 0.99789 — 08 TSO 0.99996 
20 ‘' 0.00000 0. 02374 0.29553 0. 64728 OF 80257 
0.95693 0.98924 0.99787 0.99967 0.99996 
n= 25 
ke! 1 2 5 4 5 
fo 6 fs 8 9 10. 
t 
1 ' 0.63960 0.87563 0.95906 0.98719 0.99620 
s 0.99893 0.99970 0.99991 0.99996 0.99997 
ee Bae? es. a Veta ee Cee a a a 
2 ' 0.24023 Os 78260 0.91337 0.96830 0.98931 
: 0.99667 Oe aeae OF ae 0.99991 0.99996 
5 ' 0.09031 Ue 67650 om 87087 (oe 94746 6.98052 
: 0.99343 0. oe Oe Ue a 0. ° alae 0.99994 
Seo es, | Sere” ow | eee a: ee ee ee ee | 
4 ' 0.03399 0.50116 0.81639 0. 92644 0. 97079 
: Os: es 0. suo 0. ue ole Od 0. — 
<< Eee Lee ae ee = = = pa a <<. «= to veo. ae PS 
5 0 01280 OL. 39269 oe Oe 74856 0. 90351 0.96075 
: 0. US Oe 99507 0. oes 0. oe 0.99989 
6 Oe 00483 Oe 30742 0. 68688 0. 87399 0.95043 
: Oe ” Sige 0. eee 0.99801 UF oe 0.99985 
7 Oe 00182. 0. 24078 Or 61785 0. 84046 0.93829 
: OF oe 0.99177 0. ee li}. se 0.99983 
ea ee = = ee a = = Z -— «= wo 2 Senne 
8 0% 00069. OF 18878 0. 56049 Ohe 80579 0.92397 
: 0. GOL Ox es Dak 0.99699 Oh. peel 0. eet 
Se i = = Be OS Pee aS eee See = = 
9 50. 0002 6 OP. 14822 0.50876 OF 77166 OP 90820 
: @. a 0.98864 0.99655 0.99910 0.99979 
PM of ie yd = = Se oe ee rr ee 
L0ce 0; 00010 Ge 11656 0. 46231 0.73893 0.89181 
: 0.96060 0.98704 0.99618 0.99903 0.99978 
t 
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Re 1 Z 3 4 5 
i >! 6 i 8 2 10 


i 
11 + 0.00004 0.09184 0.42071 0.70806 0.87546 
' 0. iba 0. oe 0. a 0. call 0. ee 
i = Z i ee ie = mt a s a 
12 + 0. 00001 0. 07252. 0. 38352 0. 67924 0. 85963 
' 0. pane? 0. me OF 99550 0. ae” 0393877 
Se ee re 2 é — ee ee ee 
13 + 0.00001 0. 05740 0. 35034 0. 65256 0. 84465 
' eet? 0. a 0. 0 ae 0.99890 Ue oa) 
ee es Sk ™ ss ae a = —— 
14 + 0.00000 0. 04557. Or 32081. 0.62808 0. 83078 
' 0. oS Ue 25 Ie 0: “iain 0. BOS 0.98276 
5 Se = = ae a es 
roa OF 00000 0. 03630 0. 29461. 0. 60582 0.81818 
' 0. pd 0.97786 O910'4: 2.3 0.99874 0. ood 
2 See eee cee, ome. Se = ee : 
16 «6: O. 00000 0. .02903 0.27146 Ur. 58318 2 0. 80699 
' 0. 2ohct 0. peed | 0. eyoUS Ue —— 0. ee 
= = : i 7 : 7 7 : 
3? ls «CO. 00000 OF 02333, Oe .25113 0. 56813 0. 79729 
' 0.92104 Or; tee 0529/5270 Ore sss 059277 
. 2 ae ees. ae —~ey. - = 2 eee = . 2 eee 
Be «4 0, 00000 0.01886 0.23346 OR SS283 0. 7eor 
' Oe Ss O29 7 SUS Wes 70 G.99856 2 Sooo 
ee eee a ee le |} ee 
19 + 0.00000 Oe 01537 OF 21855 0.54000 0.78252 
' Onoda 2 3 ee ool 0. pia 0.99824 0. ili 
= ORS Re ee) | ain Z owe — . = ee : nee 
20 6+ 0.00000 O20 lZos 0. 20872 0. 52969 0. 776 86 
' 0. ee 0. i 0.99184 0.99612 0.99964 
ame, ees < Ser ee — 2 Be ee ew ee See 
Zi 0. 00000 0.01055 0. 19568 0252170 0.77 aae 
t 0. hi 0. ae 0. Epa 07. ele 0. eee 
ee ae a = z 7 - = 2 — 
Eee 0% 00000 om 00898 0. 18830 Ue $1536 0. 76864 
' OF ED Oe es 0. Jone 0. sana 0. ee 
- - - t - = - - - = - = = - = = - - — = — a = 
Zs Jae 0. 00000 0. 00789 0. 18304 0. 51045 0. 76598 
' 0. 5 0. paouce 0. ee Oz ou 0. 99959 
- = - {t= = = =_ a = = - - _ — ae Ls = = 
24 + 0.00000 0. 00726 0. 17918 0. 50713 0. 16435 
' OmsGsc7 De. ee 0. ee 0. Bee 0. PES 
- - - i —_ = -_ - = = — -_ - — = — 7 a a 
25 + 0.00000 0. 00687 0. 17702 OF 50553. 0. 76367 
' 0.90564 0.96832 Oo Spe LU WS) ici 0. SSem 
t 


gs a Un SSS gg i rea fe SSS 7 SS a 


oy 





me 3 
q 
rg 6 7 
635355 0.87378 0 
. 0.99875 0.99964 
ee 23890 0.77958 0. 
0.99605 0.99877 
ee... é es 
2 ae 08947. +«:0.63210 0. 
: 0.99212 0.99737 
2 ee a 
i) 005353 0.49579 0. 
0.98732 0.99552 
a : 
2) ie 01257. 0.38683 0 
0.98196 0.99333 
a | ce ae 
6 | 0.00472 0.30143 © ©O0. 
0.97630 0.99091 
7. 0.00177 0.23488 0. 
0.97049 0.98836 
8 | 0.00067 0.18312 0. 
0.96420 0.98578 
ae a 
on ae 00025 aoeen | ala. 
0.95700 0.98313 
10 | 0.00009 0.11160 0. 
0.94886 0.98029 
11 | 0.00004 0.08726 ~=dO.. 
0.94001 0.97717 
— _— . 
i. 0. 00001. 0.06831. 0. 
0.93074 0.97376 
13 | 0.00001 0.05355 0. 
0.92134 0.97015 
14 | 0.00000 ALaloode na. 
0.91203 0.96644 
aaa ae: 
15 , 0.00000 0.03306 0. 
i 0.90301 0.96273 


.95760 0.98633 0 


QO. ee 


UG ey 


91030 Q. 96611 0. 


0. a 


86632. Q. 94367. Ue 


0299919 


0. Dae 


Ue ee 


81012 0. 92108 0. 


0. eae 


QO. mae 


"74023 QO. 89636 QO. 


0. ae 


Or: Se 


67056 0. 86460 0. 


Os amor 


0. haa 


60571. 0.82850 0. 


6.99585 


0. ak 


54670 QO. 79104 C= 


0. eS 


OQ. SEER 


(49348 0. 75399 Q. 


0. eae, 


0. oe Me 


44566 0. -71825° QO. 


Q. ha 


0. CERO! 


40277 0. “68424 Q. 


0. oo 


0. ee 


36436 Q. 65215 0. 


UE eee 


0. oe 


"32997. Q. (62203 UF: 


0. ee 


0. pe 


29922 Q. (59388 QO. 


Q. ere 


0. poe02 


27173 0. 56766 0 


0.98755 


38 


0.99638 


10 


oJ ons 
Ur 99998 
98804 
0.99996 


“97810 

Ope a 
96700. 

0. eee 
95542. 

0 eld 
94340 

0.99971 
92924 

0. oe 
“91250 

0. oe 
"89392 

99942 
874 40 
0. ee 
g5461. 

0.99925 
83503 

0.99919 
$1602 

0.99915 
79779 

0. TI 
.78049 

0.99907 





0.00000 0.02605 


0. a We NLS 


0. -00000 0. 02087 
0. poet? ee 
OF. 00000 0. 01628 
OG. sii aha 0. notes 
0. 00000 OF 01292. 
0.87246 0. aaa 
0.00000 0. 01029 


5 ee ve ee 


QO. (00823, 
a n°” UE ae 

ue (00662. 

Bes al el bas Oe Oe 
0.00009 OF 
O.635 oe 


es: 
= 


0.00000 0. 00436. 
0. $5035" 
0. 00000 Ox “00359 
On ale 


OF gee 


UF gee 
O% 00000 QO. 00299 
os —- 0. aCe 
0. 00000 QO. 00255 
Ye ane 0. gS 
0. 00000 Or 00224 
0. oes 0. aon! 
0. 00000 QO. “00206 
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